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Abstract 

The problem of finding a martingale on a manifold with a fixed ran- 
dom terminal value can be solved by considering BSDEs with a gen- 
erator with quadratic growth. We study here a generalization of these 
equations and we give uniqueness and existence results in two different 
frameworks, using differential geometry tools. Applications to PDEs 
are given, including a certain class of Dirichlet problems on manifolds. 



1 Introduction 

1.1 Martingales and BSDEs on manifolds 

Unless otherwise stated, we shall work on a fixed finite time interval [0; T]; moreover, 
{Wt)o<t<T will always denote a Brownian Motion (BM for short) in M*^^, for a 
positive integer dw Moreover, Einstein's summation convention will be used for 
repeated indices in lower and upper position. 

It is well-known that there is a deep interplay between on the one hand the 
probability theory of real martingales and Brownian motion and on the other hand 
the theory of PDEs and harmonic functions h : Mi M, defined on a manifold Mi. 
For instance, the Feynman-Kac formula gives a probabilistic interpretation for the 
solution of a PDF; besides, the Dirichlet problem for such harmonic functions h can 
be solved by considering real martingales with fixed terminal value. 

It is natural to ask whether these links can be generalized to the nonlinear 
context of manifolds, i.e. if we replace the vector space M in the definition of h 
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by a manifold M; this is what we now examine. Suppose that M is a manifold 
endowed with a connection F; then one can define the notion of F-martingale on M, 
which generalizes real local martingales (for an overview of the basic definitions and 
properties, see (TU], jS] or ^^). 

In M", the problem of finding a martingale (X()o<f<T with terminal value = U 
consists of solving the Backward Stochastic Differential Equation (BSDE for short) 



Xt+dt = Xt + ZtdWt 
Xt = U, 



where (Zt)o<t<T is a M"^'^'^ -valued progressively measurable process. With the 
connection F on M, one can define an exponential mapping exp and the equation 
under infinitesimal form (E) becomes, for martingales on M, 

(M)o( ^*+'^* = "^P^*(^*^^*) (1.1) 
I Aj^ — U 

where Zt G C{R'^^ ,Tx,M) is now a linear map. As in the linear context, studying 
martingales on M (or equivalently solving BSDE (M)o) allows to solve in a prob- 
abilistic way some nonlinear PDEs and to study harmonic mappings. Let us recall 
the definition and some properties of these mappings. 

A harmonic map H : Mi — > M between Riemannian manifolds Mi and M is a 
smooth map which is a local extremal of the energy functional 

\gTadH\\'^dvol 

where dvol is the Riemannian volume element on Mi. 

A different but equivalent point of view about these mappings is the one of a system 
of elliptic PDEs (see ji8|); let us make precise it. Consider a second-order differential 
operator C without term of order 0, defined on Mi. For h : Mi —>■ M, one can define 
by means of C the tension field of h; it is a vector field along h, i.e. 

£m(/^) : Ml ^ TM, CMih){x)en^^)M. 

Then the equation Cuih) = characterizes £-harmonic maps (see and jS], and 
probabilistic interpretations in the introductions of j2H] and [22] )■ In coordinates 
(a;*) on M and (y") on Mi, this equation can be written as the following system of 
elliptic PDEs 

V^, AM,f + F},(0)(7"^(x)D„0^D^0'= = 

with denoting the Laplace-Beltrami operator and {g°'^) the inverse metric ten- 
sor on Ml. Note that we have used the summation convention. 
With the theory of martingales on manifolds, one can solve such a system of non- 
linear elliptic PDEs; for further details, the reader is referred to and |17j . 

Now we come to the aim of this article, by enlarging the class of processes studied. 
Let (-Bl')o<t<T denote the M'^-valued diffusion which is the unique strong solution of 
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the following SDE : 

dB^, = b{By)dt + a{B^)dWt 

Bo = y, 

where a : M'^ -> M^^'^w- }, ■ _^ 

are bounded functions with bounded 

partial derivatives of order 1,2 and 3. 

On M", equation [E) is a very simple BSDE. A more general form of BSDE on 

[0;T]is _ _ _ 

studied for instance in (25 and |2S1; such an equation can be used to solve systems 
of quasilinear PDEs (see for example |25j). 

If we try to extend equation [E + D) to manifolds, by combining with equation 
(M)o we get the following equation (under infinitesimal form) 



{M + D) 



Xt+dt = exp^^ {ZtdWt + f{Bl Xt, Zt)dt) 
Xt = U. 



The aim of this work is to study existence and uniqueness of a solution to the 
generalized equation (M + D)q. 

As for martingales, which are linked to harmonic mappings h (i.e. to the above 
equation CmW = 0), this BSDE is related with a differential operator generalizing 
the tension field Cm', moreover, in some cases, the class of mappings which solve 
the new PDE can be described in terms of the local extrema of another variational 
problem. This will be discussed in Section |S] It is known that when Mi C 
and M = S"^, harmonic mappings can be used to model the state of equilibrium of 
liquid crystals (see the introduction of JH] for a brief discussion); then mappings 
associated to the new variational problem could be used to model the equilibrium 
state of a liquid crystal in an exterior field equal to the drift term / in equation 
{M + D)o. 



1.2 Setting of the problem 

In the whole paper, we will always suppose that a global system of coordinates is 
given on M . Then in these coordinates, we get from equation (M)o the following 
BSDE (see for instance JUj or the introduction of 6J) 



(M) 



dXt 
Xx - 



-- ZtdWt -\T,k{Xt){\Zt]%Zt]^)dt 



in this equation, we have used the following notations, which will be valid throughout 
the sequel : (-I-) is the usual inner product in an Euclidean space, the summation 
convention is used, and [A\^ denotes the i^^ row of any matrix A\ finally. 



/ r),(a:) \ 



1.3) 
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is a vector in M", whose components are the Christoffel symbols of the connection. 
Remark also that here Zt is a matrix in M"-^'^^. 

In this case, the classical approach of |21] to solve BSDEs with Lipschitz coeffi- 
cients fails since there is a quadratic term in Zt in the drift (the reader is referred 
to [22] or [21] for an introduction to the theory of BSDEs in Euclidean spaces). 
However, uniqueness and existence results have been obtained using differential ge- 
ometry tools, in particular by Arnaudon (PP), Darling ([!]), Emery(|lOj). Kendall 
([T3]). Picard ([2Z| and [21]) or Thalmaier ([^21 and [SI]); note also the results of 
Estrade and Pontier (pi]) concerning some classes of Lie groups. Independently of 
geometric tools, a lot of works have tried to weaken the Lipschitz assumption : in 
the one-dimensional case, they include [18] (in dimension one, her results are more 
general than the ones of this paper because she deals with generators with quadratic 
growth), [20] or ^2]; in higher dimensions, we refer the reader for instance to |3], 
PU] (who studies a Ricatti-type BSDE) and [22]. To the best of our knowledge, 
there is no paper that would include our results in dimensions greater than one. 

Now in our global chart, the equation (M + D)o becomes 

(M + D) ( = + {-\V^,{Xt){\ZtY\[Zt]^) + /(i?f ,X„ Zt)) dt 

(the same notation will be used to denote the TM-valued function / and its image 
in local coordinates). The process X will take its values in a compact set, and a 
solution of equation {M + D)o will be a pair of processes (X, Z) in M x (M'^w-- ^tM) 

such that X is continuous and E ^ \\Zt\\1dt^ < oo (|| ■ ||^ is a Riemannian norm; 

see below). If we consider a global system of coordinates on an open set O of R", it 
corresponds to processes (X, Z) in O x W"^^ , such that X is in a compact set and 

E ^ WZtW^dt^ < oo, solving equation (M + D). 

Two different cases will be considered here : firstly when the drift / does not 
depend on and secondly the case of a general / in nonpositive curvatures. In the 
two Riemannian structure is fixed on M; in the former case, the connection 

may be independent of this Riemannian structure, while in the latter, only the Levi- 
Civita connection associated will be used. Note that the case for a general / with 
K > Q involves more technical calculations; it will appear elsewhere. 

We first give in Section |21 mild generalizations of well-known results, concern- 
ing the geometry of the manifold and a characterization of the solutions in the 
z-independent case by means of convex functions. In Section |21 we study the unique- 
ness problem. It is solved by generalizing to our context two methods : on the one 
hand, Emery's idea, used in pQ]j and [15 ; on the other hand, the work of Picard 



We obtain Theorems 13 . 3 . 21 and . 4 . fil 
Section m is devoted to proving the existence of a solution of equation (M + D). 
The main arguments are to exhibit a solution for " simple" terminal values (based on 
a strong bound on the process (ZJ in Subsection 14. 3|) and to solve the equation for 
any terminal value using approximation procedures (Subsections 14.11 and I4.6|l . We 
need for the proof an additional (and necessary in fact) condition on the drift / : 
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it is supposed to point outward on the boundary of the set on which we work. We 
give in Subsection 11.41 the main result (Theorem 11.4.11) which sums up the results 
obtained. In Sectional we extend the results to random time intervals [0; r], where 
T is successively a bounded stopping time (Theorem 15.3. ip and a stopping time 
verifying an exponential integrability condition f Theorem I5.3.2j) : then to conclude 
this paper, we give some generations and applications to the theory of PDEs, as well 
as the variational problem related to equation (M + D)q. 

The uniqueness part, as well as the applications to PDEs, are mainly adaptations 
of procedures already used; on the contrary, the approach for the existence seems 
to be novel. 



1.3 Notations and hypothesis 



In all the article, we suppose that a filtered probability space {fl,J-',P, (J-'t)o<t<T) 
is given on which {Wt)t denotes a rfviz-dimensional BM. Moreover, we always deal 
with a complete Riemannian manifold M of dimension n, endowed with a linear 
symmetric (i.e. torsion-free) connection whose Christoffel symbols F*^ are smooth; 
the connection does not depend a priori on the Riemannian structure. 

On M, 6 denotes the Riemannian distance; \u\r is the Riemannian norm for a 
tangent vector u and \u'\ the Euclidean norm for a vector u' in M". If /i is a smooth 
real function defined on M and Mi,M2 are tangent vectors at x, the differential of h 
is denoted by Dh{x) < u > or h\x) < u >; the Hessian Hess h{x) is a bilinear form 
the value of which is denoted by Hess h{x) < ui,U2 >■ 

For /? G N*, we say that a function is on a closed set F if it is on an 
open set containing F. For a matrix z with n rows and k columns, *z denotes its 
transpose. 



J2i=i \ [^^Y\r "where the columns of 



(Tr is the trace of a square matrix) and \\z\\r = 
z are considered as tangent vectors. The notation \E'(x, x') ~ 6{x,x') 
there is a constant c > such that 



means that 



Vx, x', - 6(x, x'Y < \l/(x, x') < c 6(x, x'Y. 
c 

Finally, recall that a real function x defined on M is said to be convex if for any 
M-valued geodesic 7, x ° 7 is convex in the usual sense (if x is smooth, this is 
equivalent to require that Hess x be nonnegative). 

Before the general framework, let us give some additional notations which are 
specific to the Levi-Civita connection. In this case, we always assume that the 
injectivity radius i? of M is positive and that its sectional curvatures are bounded 
above; we let K be the smallest nonnegative number dominating all the sectional 
curvatures. For the distance function 6, if x = {x, x') is a point and m, u are tangent 
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vectors at x', we consider the partial derivatives (when they are defined) 

62{x) < u >= 6'{x) < (0,m) > 

and 

Hess22 5(5) < u,u >= Hess S{x) < {0,u), {0,u) > . 

Now let us recall from ^3] the definition of a regular geodesic ball. A closed geodesic 
ball B of radius p and center p is said to be regular if 

(i) pVk<^ 

(ii) the cut locus of p does not meet B. 

For an introductory course in Riemannian geometry, the reader is referred to and 
for further facts about curvature, to P^ . 

Throughout this article, we consider an open set a; 7^ relatively compact in 
an open subset O of M, such that there is a unique 0-valued geodesic between any 
two points of O; O is also supposed to be relatively compact in a local chart, so 
that it provides a system of coordinates (global on tJ); it will be as well considered 
as a subset of M". We suppose that there exists a nonnegative, smooth and convex 
function \l/ on the product tJ x cJ (i.e. convex on an open set containing this set) 
which vanishes only on the diagonal A = {{x,x)/x &uj} {uJ is said to have F-convex 
geometry); besides, we suppose that ^ 6^ for ap > 2 (note that since is smooth, 
p is an even integer). In fact, we take for ZJ a sub level set of a smooth convex function 
X defined on O : {x < c} (note that this hypothesis and the existence of \Ef guarantee 
the existence and uniqueness of a cJ- valued geodesic between any two points of ZJ). 
Emery has shown (see Lemma (4.59) of 10 ) that, in the case of a general connection 
F, any point of M possesses a neighbourhood with F-convex geometry; when the 
Levi-Civita connection is used, this is true for a regular geodesic ball (see piHj). 

Finally we always assume two hypothesis on / : 

3L > 0, \/b,b' G M^V(x,z) G O X /:(M'^^,T,M),V(x',;z') G O x /:(R'^^, T,,M), 



\fib,x,z)-fib',x',z') 



<Li{\b-b'\ + 6{x,x')){l+\\z\\r+\\z'\\r) 



+ 



\z — z 



;i-4) 



and 



3L2 > 0,3x0 G 0,V6 G R"', 1/(6, xo, 0)1, < L2. 



;i.5) 



The first one is a "geometrical" Lipschitz condition on /. This special form is needed 
to get an expression which is invariant under changes of coordinates. We will see 
that later, in ()2.5j) . Remark that, in the 2;-independent case, it just means that / is 
Lipschitz with respect to the first two variables; otherwise. 
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denotes the Riemannian parallel transport along the unique geodesic between x and 
x' . The second one means that / is bounded with respect to the first argument. 
Remark that these conditions also imply the boundedness of / if it does not depend 
on the variable z. 

Note to end this part that the same letter C will often stand for different constant 
numbers. 

1.4 The main result 

Before achieving calculations, we give the main theorem of the article. Let us first 
introduce a technical but natural hypothesis, which we will make explicit in Subsec- 
tion SH: 

{H) f is pointing outward on the boundary of ZJ. 
Then we can state : 

Theorem 1.4.1 We consider the BSDE (M + D) with terminal random variable 
U & uj = {x < c}, where oJ satisfies the above conditions. If f verifies conditions 
(jl.4|) . (|1.5|) and {H), and if x is strictly convex (i.e. Hess x is positive definite), 
then 

(i) If f does not depend on z, the BSDE has a unique solution {Xt, Zt)o<t<T such 
that X remains in uJ. 

(a) If M is a Cartan-Hadamard manifold and the Levi-Civita connection is used, 
then the BSDE has yet a unique solution {Xt, Zt)o<t<T with X in uJ. 

In particular, if the Levi-Civita connection is used, a "good" example of domain on 
which existence and uniqueness hold is a regular geodesic ball. 

In Section we will extend this theorem to random time intervals [0; r] (in- 
stead of [0;T]), for stopping times r which are bounded, or verify the exponential 
integrability condition : 

3p > : E(e^O < oo. (1.6) 

In the former case. Theorem 11.4.11 goes the same, while in the latter the constants 
L and L2 in ()1.4|) and ()1.5|) are furthermore required to be small with respect to the 
constant p in p.6|) . 

Acknowledgements : The author would like to thank his supervisor Jean Picard 
for his help and his relevant advice, and the referees for their suggestions to improve 
a first version. 

2 Preliminary results 

We first recall elementary results about Ito's formula and parallel transport. Then 
we give some geometrical estimates for the distance function on M x M and char- 
acterize solutions of the equation (M -|- D) using convex functions, but only when 
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the drift / does not depend on z. As underlined in the introduction, these results 
are just mild generalizations of well-known results of ^TU] and [2H]. 

In this section, the covariant derivative of a vector field Zt along a curve 7i will 
be denoted V^^-Zf. 

2.1 Ito's formula on manifolds 

Consider two solutions and (X^, Z^) of equation (M + D) with terminal 

values f/^ and f/^, such that and remain in O . Let 

X=(X\X2) and Z = 

then Ito's formula with the function \E' is written 




^(Xi)-^(Xo) = y D^(x,) (z.rfiy, 

^Jo ^ V fiBy^xiz',) - '^T,,{x!){[z^nz^Y) J 

Tr (%D2^(X,)Z,) 



+ 2 I (E*[*^^]'HessM/(X,)[*Z,]Mrfs 



,«=i 

(remember notation (jl.Hp and that ['^Zsf denotes the i^^ column of the matrix Z; it 
is a vector in R^"). Moreover, for a smooth function h on O and a solution (X, Z) 
oi {M + D), we get a similar formula, replacing X by X and Z hj Z. 



2.2 Two inequalities 

We first give an equivalence result between the Euclidean and Riemannian norms; 
it follows easily from the relative compactness of O, considered as a subset of M" (in 
particular, this means that we can identify each tangent space with M"). 

Lemma 2.2.1 There is a c> such that for any (x, z) e O x T^M = O x M", 

1 

-\z\ < \z\r < C\z\. 

c 

Lemma [2 . 2 . 1 1 will often be useful in the sequel. 
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Proposition 2.2.2 The Levi-Civita connection is used. There is a C > such 
that for every {x, x') E O x O and {z, z') G Tr^M x T^'M, we have 



and 





x' 








\\z-z' 


< 


C{\z- 




X 


r 








x< 


\z 


-z'\<C\ 


\\z-z' 






X 



+ 6{x,x'){\z\r + \z'\r) 



(2.2) 
(2.3) 



Remark 1. In fact, by Lemma f2.2.H we can use any of the two norms (except 
for \z — z'\ which is necessarily the Euclidean norm). 
Proof. It is sufficient to prove 



y{x,z),{x',z') eOx 
Indeed, this implies 



\z — z 



< C5{x,x')\z\ 



(2.4) 



\z — z 



<Cii\z- z'\ + 



\z — z 



< C{\z- z'\+6{x,x')\z\) 





x' 




x' 




x' 




\z - z'\ < 


\\z — z 

X 


+ 


\\z-z' 

X 


< C [5{x,x')\z\ + 


\\z-z' 

X 


) 



and 



So let us prove ()2.4|) : let 7 be the geodesic such that 7(0) = x and 7(1) = x', 
and z{t) the parallel transport of z along 7 : 

VtG[0;l], z{t)=\z. 

X 

In local coordinates, the equation V^t-2(t) = gives for every k 

z\t)+T%{^,nz\t) = ^. 

Moreover, \'^t\r = 5{x,x') so |7f| < C5{x,x') and 



\z — z 



E 

k 



Z^{t)dt 



T'^,{^,)^lz\t)dt 



< Ci6{x,x'fy2 f \z\t)\^dt 
, Jo 



< C6{x,x'y\z^^ 
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The last inequality comes from the equivalence on O of the Riemannian and Eu- 
clidean norms, and the fact that is constant by definition of z{t). The proof 
is completed. □ 

As a consequence, on the relatively compact set O C M", ()1.4|) becomes 
3L' > 0, V6, b' e W^, V(x, z), (x', z')eOx M""^*^, 



1/(6, x, z) - f{b', x\ z')\ < L\{\b - b'\ + \x- x'|)(l + ll^ll + ll^ll) + \\z- z'W). (2.5) 

Remark 2. In particular, ()2.5p is verified for a drift / that is Lipschitz in 
(6, X, z) in O (the Lipschitz property of / is not necessarily preserved by a change 
of coordinates, but (j2.5|) is). 

Remark 3. If / does not depend on it just means that / is Lipschitz in (6, x). 



2.3 Estimates of the derivatives of the distance 



This paragraph is based on Section 1 of j2Hl- The connection used is Levi-Civita's 



one. 



The geodesic distance (x, x') i— > (5(x, x') is defined on M x M and is smooth 
except on the cut locus and the diagonal {x = x'}. We want to estimate its first and 
second derivatives when M x M is endowed with the product Riemannian metric. 
If X = (x, x') is a point which is not in the cut locus or the diagonal, there exists 
a unique minimizing geodesic 7(t), < t < 1, from x to x'. If Ut is a vector of 
T^(^t)M, we can decompose Ut as Vt + Wt, where Vt is the orthogonal projection of Ut 
on 7(t); the vectors Vt and Wt are respectively called the tangential and orthogonal 
components of Ut- If m = (mq, Mi) is a vector of T^{M x M), (fo, f i) and {wq, Wi) are 
also called its tangential and orthogonal components. 

Lemma 2.3.1 Let x he a point of M x M which is not in the cut locus or the 
diagonal. Let u be a vector of Tx{M x M) and let v and w be its tangential and 
orthogonal components. Then 



\S'{x) <u> 



\Vo-Vi 



(2.6) 



if moreover K = (i.e. the sectional curvatures are nonpositive), then 

2 



Hess (5(x) < u,u>> 



S{x) 



\Wo-Wi 



(2.7) 



Proof. Let Jv{t) (resp. Jw{t)) be the tangential (resp. normal) Jacobi field along 
7(t) satisfying J„(0) = vq and Ji,(l) = vi (resp. Jw{0) = wq and Juj(l) = wi). From 
(1.1.5) and (1.1.7) of |2H1, we have 



6'(x) <u>-- 



(7(t)|V^(f)J.(t)) 

m)\r 



{2.i 
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and 



1 

Hess 6(x) <u,u >> I 



V^it)Mt)\idt-K6{x) / \Mt)\idt. (2.9) 



It is easy to see with the Jacobi equation that Jv{t) = {At + B)'y{t). Then 
V-y(t)Jv{t) = A'jit) and the hmit conditions J^(0) = vq = 07(0) and Jv{l) = vi = 
/37(1) imply A = jd — a. Hence ()2.8p gives 



\5'{x) <u>\ = \A\- |7(t)|r- = 1/3 - a\5{x) 



\Vq-Vi 



This is 

Moreover, let us write Jwif) = X]?, '"«('^)-^«(^) where {Ei(t)}i is a parallel or- 
thonormal frame along 7. Its covariant derivative is V.y(t)Jw(t) = J2i'^ii't)Ei(t) and 

\\MO) = EiUMH^)- Then 



2 2 

^(0) - JM) 



Now this inequality together with p.9|l and the nonpositivity of the sectional cur- 
vatures give ()2.7j) . □ 

Then we have the following estimate : 
Proposition 2.3.2 If K = and x is not in the cut locus, then 

,1 



Hess {^S'^){x) <u,u >> 



I Mo - Ml 



:2.10) 



Proof. If X = {x, x') is not on the diagonal, we recall the classical formula 
Hess <u,u >= 6{x) ■ Hess 6{x) <u,u> +{S'{x) < u >)^. 

This formula and estimates ()2.6|) and ()2.7p imply the proposition for x ^ x' since 

x' x' 

the two terms ||mo — Vi and \\Wq — Wi are orthogonal for the Riemannian norm. 

X X 

The follows by continuity since 5^ is smooth on a neighbourhood of the 

diagonal. □ 
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2.4 A characterization of the solutions of equation (M + D) 
when the drift / does not depend on z 

We give here a generalization of a well-known result (see (4.41) (ii) in ^Uj) which 
roughly says that a continuous M-valued process {Yt) is a F-martingale if and only 
if its image under convex functions is a real local submartingale. In this paragraph, 
the filtration used is the natural one of {Wt)t- 

Proposition 2.4.1 Suppose that the drift f does not depend on z. Then every 
point p of M has an open neighbourhood Op, included in a local chart, with the 
following property : 

A pair of processes {X,Z) (with X continuous, adapted and Op-valued) is a 
solution of (M + D) iff for every convex function ^ : Op — » M, i{Xt) — /q D^{Xs) ■ 
f{B^,Xs)ds is a local submartingale. 

The proof given here is just an adaptation of Emery's one; first we recall Lemma 
(4.40) of P]. 

Lemma 2.4.2 On M, let C, be a smooth function. Every point of M has an open 
neighbourhood Op depending on ^ with the following property : For every e > and 
a G Op, there is a convex function h'^ : Op ^ W such that {a,x) ^ h'^{x) is smooth 
on Op X Op, h'^{a) = 0, Dh'^{a) = -D^(a) and Hess h'^{a) = eg{a), where g represents 
the metric. 

Now we complete the proof of Proposition ()2.4.1|) : 
The "only if part is just a consequence of Ito's formula (similar to ()2.H) ) applied 
to ^{Xf) : as Hess ^ is nonnegative by convexity, C,{Xt) — /g DC^{Xs) ■ f^B^, Xs)ds is 
indeed a local submartingale. 

For the "if part, notice first that around p there is a system (x*) of local coor- 
dinates that are convex functions (if (y*) are any local coordinates with y\p) = 0, 
then take = + cJ^jiV'')'^ for c > large enough). Choose Op relatively compact 
in the domain of such a local chart and in an open set on which Lemma ()2.4.2j) holds 
for ^ = X* and ^ = — x* (denote by 5* this set of 2n functions). 

For a continuous adapted Op-valued process X, suppose that ho X — f Dh{X) ■ 
f{By,X)dt is a local submartingale for every convex h on Op. Taking first for h 
the global (on Op) coordinates (x*) shows that each x* oX is a real semimartingale, 
so X is a semimartingale. For every fixed ^ G S", it is sufficient to prove that 
^oX-|/(Hess ^)ijd < X\Xi > - J Di^{X) ■ f^{By , X)dt is a local submartingale; 
for then replacing ^ by — ^ shows that (remember that DijX^ = and D^x^ = 1 if 
i = k and otherwise) 

x^ o X + ^ j T'yjix)d < x\ X' >- j f\By, X)dt 

is a local martingale for each k. But the theorem of representation of local mar- 
tingales in Brownian filtrations allows to write this local martingale explicitly as 
/ ZtdWt, thus d < X\X^ >t= {[Ztf\[Zty)dt and (X, Z) solves equation (M + D). 
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By the choice of Op, given any £ > we are provided with functions h'^ associated 
to ^ as in Lemma ()2.4.2|) . Call a the p^^ dyadic subdivision of the time axis (cr = 
{tfc = |r : k,p e N}) and let for t G [tk] p{t) = tk and 

St = E f - (^") ■ fi^u, Xu)du) 

As each /i" is convex, he (Xt) — j^^ Dhe {Xy) ■ fiB^, Xu)du is a submartingale in 
the interval [t^S^+i], and S'' is a continuous submartingale. Using the coordinates 
(x*), write 

dS^ = D,he''{Xt)dXi + '^■(X^d < X\X^ >t ~D,hf''{Xt)f{BlXt)dt. 

So, if X* is decomposed into A^* + (i.e. local martingale + bounded variation 
part), 

dS^ - (^D,hf''{Xt)dAi + ^D,jh^'^{Xt)d < X\X^ >i ~D,h'^''^{X,)f\BlXt)d?j 
is a local martingale; hence the process 

B^, = f D,h^'''\Xs)dAl + l f D,,hf'''\X.,)d<X\X^ 
Jo ^ Jo 

- [ D,hf'''\X,)-riBy,Xs)ds 
Jo 

is increasing. The estimates (using relative compactness of Op) 

\Dih^^{v) - Dihl{v)\ < C\u-v\ 
\D,,h^{v) - D,,h^,{v)\ < C\u-v\ 



and the convergence of A^p(t) to Xt when p goes to infinity (i.e. when a becomes 
finer), yield a dominated convergence of Dike ''^■^\Xs) to Dihf''{Xs) = Di^{Xs) 
(since rf/i^(a) = di{a)) and of Aj/^f (X,) to rf^.(X,)Dfc^(X,) + egij{X,) (since 
Hess /i"(a) = £:(?(a)). Hence -8°" has a limit, equal to 

j D,i{X)dA' + ]^jT\^{X)Dki{X)d<X\X^> 

+ \^ j 9^j{x)d < x\x^ >- j AeW ■ f\By,x)dt, 



13 



that is an increasing process too. Letting now e tend to zero, 
J = j D,i{X)dA^ + ]^ j V\^{X)Dui{X)d<X\X=>- j D,i{X)-r{By,X)dt 
is also increasing, and (remember that Dij^ = 0) 

e o X - e o Xo - i y (Hess 0.,^^ < X\ X^ > 

- J Di^{X) ■ f\B\ X)dt = j Dii{X)dN' + J 
is a local submartingale, as was to be proved. □ 

3 The uniqueness property 

In the first paragraph, we set the problem and exhibit the sum (j3.2|) . whose non- 
negativity suffices to have uniqueness. Then, we give a useful estimate and derive 
the result in the two cases considered in this paper. The calculus is rather longer if 
the drift / depends on for we have to prove exponential integrability. 

3.1 The general method 

Consider two solutions (X^, Zt)t and (X^, Z'^t of {M + D) such that X and X' remain 
in uJ and Xt = Yt = U (we will often write "cJ-valued solutions of (M + -D)"). Let 



X, = (X„X;) andZ, = (^^^^ 



In the martingale case (/ = 0), Ito's formula ()2.Hl and the convexity of \Ef ensure 
that the process (\E'(X())t is a submartingale. The other properties of \E' (see the 
introduction) then imply that X remains in the diagonal A, therefore the uniqueness 
required. This is Emery's method (see Corollary (4.61) in [TU]). 

For our purpose (i.e. / does not vanish identically), we want to keep the sub- 
martingale property and therefore control the integral involving / in p.lj) . The idea 
is to study, rather than (\I'(Xj))f, the new process (exp(y4t)\I'(Xf))t where 

At = Xt + iJ, {\\Zs\\r + \\Z'J\r)ds, 
Jo 

for appropriate nonnegative constants A and fj,. Apply Ito's formula to obtain 
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Jo Jo 



fiBy,xi,zi: 

+ / e^^^(X,)(A + Mll^s||r+||^:i|r))c?s. (3.1) 



It is clear that the submartingale property will be preserved if we show the 
nonnegativity of the sum 



i=l ^ ^ ti t 

+ {X + fi{\\Zt\\r+\\Z',\\r))^{Xt). (3.2) 

The remainder of the uniqueness part is mainly devoted to proving this result. 
3.2 An upper bound 

The first step towards the nonnegativity of the sum ()3.2|) is to give a bound on the 
term involving /. This is the purpose of Lemma [3.2. II below. 

Let {x, x') be a point in aJ x cJ and b G M'^. For notational convenience, we 
keep the same notation aJ x cJ for the image of this compact set in local coordinates 
considered below, and we write / for f{b, x, z) and /' for /(6, x', z') (note that the h 
is the same in / and /'). Take a local chart (0,0) in which (x, x') has coordinates 
(£, x'); if ((9i, . . . , d2n) denotes the natural dual basis of these coordinates, then 



^ ^ 1=1 



(3.3) 



Now V = {vi, . . . ,V2n) = {x — x',x') are new coordinates in which the diagonal A 
is represented by the equation {vi = ■ ■ ■ = f„ = 0}; moreover, if {ipi, . . . , ip2n) is 
the natural dual basis of the f-coordinates, then ipi = di and ipi+n = di + di+n for 
z = 1, . . . , ra; thus ()3.3|1 becomes 



i=l 

and 



^ / 7 — 1 



15 



Finally, we call piV) the projection of a vector V onto A : ii V = {vi, . . . , V2n), then 

P{y) = (0,...,0,t;„+i,...,W2n). 

Lemma 3.2.1 Suppose that '^{x,x') ~ 5{x,x'y on uj x U where u is an even 
positive integer (since is smooth). Then there is C > such that, for all b in M.'^, 
X, x' in To and z, z' in W'^'^'^ 



f 
f 



<C6{x,xr~HSi^,x')\f\ + \f-f'\)- (3.5) 



Proof. First remark that on u x uj, we have 

\V -p{V)\^ 6{x,x'). (3.6) 

Indeed, let {x,x') G M x M; in f -coordinates, it is represented hj V = {x — x',x'). 
Then |^ — p(V^)| = \{x — x',0)\, the Euclidean distance between x and x'; as x and x' 
belong to the compact set aJ, the Euclidean and Riemannian distances are equivalent 
and ()3.6p follows. 

Write Di for Since ^{x,x') ~ S{x,x'y, a Taylor expansion gives 
^(^) = ^ E D,^...^.^{p{V)){V~piV)),, . . . (V-piV)),^ + 0{\V-piV)r'). 

l<ii,...,iv<n 

Using another Taylor expansion with Dj\E' gives 
|A^(^)| < C\V - p{V)\''-'^ iii<n 

|Dj\l/(V^)| < C\V — pCV)]" if i > n (i.e. in the direction of the diagonal) 
with a uniform C on the compact set cJ xtJ. As a consequence of ()3.4j) . we obtain 



/ 



<C\V-piV)r'{\V-p{V)\-\f'\ + \f-f'\)- 



Using once again ()3.6|) completes the proof. □ 

We now give two functions which verify the hypothesis of Lemma ()3.2.ip . In 
these two examples, the Levi-Civita connection is used. 

Example 1. Consider the mapping \E'(x,a;') = 6^{x,x'); it is smooth, not convex 
in general, but this is true if the sectional curvatures are bounded above by 0. 

Example 2. Take for uJ a regular geodesic ball B centered at o G M, with the 
sectional curvatures bounded from above by a constant K > 0; then we can find a 
nonnegative convex function on B x B which vanishes only on the diagonal (see 

m) ■■ 

^ ^ f 1-cos{VK5{x,x')) 

ycos(v^5(x, o)) cos(v^5(a;', o)) - /i2 

where /i > is small and p > 2 is an integer large enough (so that ^ is smooth). 
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3.3 The case / independent of z 

In all this paragraph, / does not depend on i.e. f{b,x,z) = f{b,x). 

Proposition 3.3.1 If the drift f doesn't depend on z, then the process [e'^^^{Xt))t 
is a suhmartingale for A positive large enough; this implies that two uJ-valued con- 
tinuous semimartingales (Xt) and (Xj) verifying the same equation (M + D), with 
the same terminal value, are indistinguishable. 

Proof. For x,x' G aJ, \f{b,x) — f{b,x')\ < L'6{x,x'); thus using ()3.5|) and the 
boundedness of / 

Then for A > C, /i = 0, the sum (|3.2p is nonnegative (note that the convexity 
of gives the nonnegativity of the term involving Hess \&). Moreover, the local 
martingale in equation 1)3.11) is in fact a martingale since D^^ is bounded on tJ x U, 
so the process (e'^*\E'(Xt))t is indeed a submartingale. As it is nonnegative and 
has terminal value 0, it vanishes identically; so \E'(X() = for all t. Finally, the 
definition of \Ef leads to Xt = X[ for all t and the proof is completed since we 
consider continuous processes. □ 

Remark. Of course, X = X' implies that for any t, Zt = Z[ a.s. 

As an immediate corollary, we are now able to give the uniqueness property. 

Theorem 3.3.2 Suppose that uJ and \E' verify the properties of the introduction 
(see paraaraph M.'Jji and moreover that the drift f depends only on {b,x) and verifies 
()1.4|) and ()1.5|) . Then, for a given terminal value U in uJ, there is at most one 
UJ-valued solution to the equation (M + D). 

Example. Suppose that the Levi-Civita connection is used. Then Theorem 
13.3.21 implies uniqueness on any compact set of a Cartan-Hadamard manifold, and 
on any regular geodesic ball (with K > Q)] indeed it suffices to consider respectively 
the functions 5"^ and ^ defined after Lemma [3.2. II 

3.4 The general case in nonpositive curvatures 

In this subsection, the drift / depends also on 2;, \I' = and the connection used 
is Levi-Civita's one; moreover M is supposed to be a Cartan-Hadamard manifold 
(i.e. simply connected with nonpositive sectional curvatures); remark then that any 
closed geodesic ball is regular. By achieving explicit calculations we are going to 
derive the uniqueness property for any compact set. 

The problem is to show that the process (exp(74t)\E'(Xj))t is a submartingale. But 
to define such a process, we need to consider solutions in some class which we now 
define. 



< C5{x,x'y < C'${x,x'). 
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Definition 3.4.1 If a is a positive constant, let (Sa) be the set of uo-valued solutions 
of {M + D) satisfying 



T 

2, 







Eexp a / \\Zs\\;ds < oo. (3.7) 



Actually, we now verify that for a small, {Sa) contains any solution of equation 
(M + D). The first step is the following lemma, which generalizes Proposition 2.1.2 
of 



Lemma 3.4.2 Suppose that we are given a positive constant a and a function 
(ponuj satisfying Cmin < < Cmax for some positive Cmm md Cmax- Suppose 
moreover that Hess (p + 2a(j) < on To; this means that 

Hess <u,u> +2a(j){x)\u\l < 0. (3.8) 

Then, for every e > 0, any U-valued solution of (M + D) belongs to (Sa-e)- 

Proof. Define 

St = (j){Xt)exp(a [ WZsWlds--^ [ W^X,) ■ f {8^, X^, Z,)\ds 

V Jo '^min Jo 

Denote by e{t) the exponential term above. It follows from Ito's formula that 
dSt = e{t) ( D(l){Xt){ZtdWt) 

dw 



+ J]*[*Zt]'Hess 0(Xj[%]* + a0(Xi)||Z,||2j dt 

+ {-'^WiXt) ■ f{BlXt,Z,)\+D<p{X,) ■ f{BlX„Z,)) dt 

Thus condition ()3.8|) ensures that St is a local supermartingale; since it is nonneg- 
ative, we get E^t < ESq. By using the lower and upper bounds on 0, we deduce 
that 

Eexp/'a \\Z,\lds--^ f \D<\>{X,) ■ f{Bl,X,,Z,)\d^ ^C^_ ^gg^ 

V Jo ^min Jo J ^min 

By conditions ()1.5|) and ()2.5|) . there is a C > such that 

-^W{Xs)-f{By,Xs,z,)\ < ^(i + iiz.ii,) 

< e\\Zs\\l + Ce; 

in particular, we have (|3.7|) for a — e. D 
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Lemma 3.4.3 Let B{o,p) be a (regular) geodesic ball of center a and radius p. 
Then for a > small enough, every B{o, p)-valued solution of (M + D) belongs to 
the set {£a) defined previously. 

Proof. Let = cos ^^(5(o, x)j for x in the geodesic ball B{o,p). We want 

to prove that this function verifies the hypothesis of Lemma ()3.4.2|) . Firstly, is 
obviously a smooth function and there is a c > such that 



Vx e B{o,p), c < (j){x) < 1. 

Moreover, we have 

Hess (j)(x) < u,u >= — cos ( — 6(o, x)] ( — ) (SUo, x) < u >)^ 

\3p J \3pJ 



IT 



( 



— sin ( — (5(o, x) ( — Hess22 5(o, x) < m,m > . 
V3p ) \3pJ 

If we get back to the notations of Subsection 12.31 we have, as a consequence of ()2.6|) 
and ()2.7p . the estimates (for x 7^ o) : 

1^2(0, X) < U > \ = \v\r 

and 



\w\ 



}less22 S(o,x) < u,u >> ^, ,. 

0(0, x) 

Then 



Ress<P{x)<u,u> < -0(x)( -J \v\l- ^^^^^^^ i^-j HI 



3p 



2 



So there is a a > (depending on the radius p) such that ()3.8j) holds. It suffices to 
apply Lemma l!l4.2l to conclude. □ 

Obviously, this is also true for any compact set aJ. We can easily derive, from 
()3.9|) and the preceding lemma, the following result, which will be useful in later 
calculations. 

Corollary 3.4.4 There is an a > and a finite positive constant Cu (both de- 
pending only on uJ and the constants L and L2 in ()1.4|) and ()1.5p ) such that for any 
W-valued solution (X, Z) of equation (M + D), 

Eexp (a [ \\Zs\\lds] < C^. 
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We can now state the uniqueness result f Theorem I3.4.6|) . 



Lemma 3.4.5 For two solutions (X, Z), {X', Z') of {M + D), verifying E(At) < 
oo for every fi > 0, the expression ()3.2|) is nonnegative for A and ^ large enough. 



Proof. Using ()2.1()j) we have for 



z 
z' 



d-w 



^*[*5]^Hess ^{x)f~zX > J2 



i=l 







2 


x' 










II I'zT - 


- h'\ 




\\z-z' 






r 


X 



(3.10) 



Moreover, we also have using p.4j) and ()1.5|1 . together with ()3.5|1 (or with ()2.2|1 
and (I21I)) 



f{b,x,z) 
f{b,x',z') 



< Ci6{x, X') ( 6{X, X'){1 + \\z\\r + Wz'Wr) + 

1 



\z — z 



< C(5'(X,X')(1 + ll^llr + Ik'llr) + 



\z — z 



(3.11) 



Then ()3.2|) is greater than the following sum 

-C6\Xt,X',){l + \\Zt\\r + \\Z',\\r) + ^6\Xt,X',){X + fl{\\Zt\\r + \\Z',\\r)). 

Taking A and fi greater than 2C makes obviously this sum (and the expression 
(I3.2|l ) nonnegative. □ 



Theorem 3.4.6 Suppose that M is a Cartan-Hadamard manifold and that the 
drift f verifies condition ()1.4j) and ()1.5|) . Then, for a given terminal value U in 
the compact ZJ, there is at most one TD-valued solution to the equation (M + D) 
(i.e. for two solutions (X, Z) and (X', Z'), the (continuous) processes X and X' are 
indistinguishable) . 

Proof. Every compact uJ is included in a closed geodesic ball. Hence any cJ- valued 
solution of (M + D) is in (Sa) for a a > from Lemma (3.4.31 this easily gives the 
integrability of exp(y4j') for every /i > 0. Consequently, it suffices to apply Lemma 
13.4.51 and conclude as in the proof of Proposition 13.3.11 □ 

Remark. Note that the only hypothesis required for uniqueness in this case is 
the compactness of ZJ. 
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4 Existence results 



In this section we are given an cJ- valued random variable U and we want to construct 
a pair of processes (X, Z), satisfying equation (M + D), with X in cJ and terminal 
value U . We limit ourselves to the case of a Wiener probability space and we recall 
that on tJ, if (X, Z) and (X', Z') are two solutions of the equation (M + D) and 
X = (X, X'), then for A > and /i > large enough the following processes 

(e^'^iXt)) iif = f{By,X) 

V /tG[0;T] 

and, when M is Cartan-Hadamard, 

V /tG[0;T] 

are nonnegative submartingales. 

The strategy of the proof can be described as follows : 

1. Simplify the problem by considering only terminal values which can be expressed 
as functions of the diffusion at time T, i.e. U = F{Bj) f Subsection 14. This 
step needs to pass through the limit in equation (M + D)] when / is independent 
of z, it is a corollary of a well-known result, but in the other case, more technical 
calculations, using the uniqueness part, are involved. 

2. Solve a Pardoux-Peng BSDE with parameter to construct a pair of processes in 
X W"^"^ which is close to being a solution of (M + D) with Xt = U (Subsection 

US). 

3. Show that under an additional condition on / the solution of the preceding BSDE 
is a solution of the BSDE (M + on a small time interval (Subsections 14.31 and 
I4.4|l . Note that the main argument in the existence proof is certainly Proposition 
14.3.21 where we give an a.s. upper bound on the process {Zt). 

4. Use the convex function \E' to show that we have a solution of (M + D) on the 
whole time interval [0; T] (Subsection 14. 5|) . 

In fact, for technical reasons we suppose in the two last steps that / is sufficiently 
regular; then the proof of the existence is completed with the last subsection : 

5. Solve BSDE (M + D) for general / using classical approximation methods (Sub- 
section H^I)- 

Note that we usually work within local coordinates in M", i.e. we consider that 
cJ C O C M". 

4.1 Reduction of the problem 

Let C^(R'^,a7) denote cU- valued functions on which are constant off a compact 
set. In this paragraph, it is shown that it suffices to check the existence result for 
U = F{Bl) with F e C^{W^,uo). 

The space of all functionals {G{Wt^,Wt^, . . . , W^), < h < . . . < tq < T; G e 
C^{W^'^,uj)} is dense in L^{Tt;uj) (remark that L'^{J^t;uJ) = L°°(J^T;aJ) since aJ is 
compact) endowed with the distance 
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The same is true by replacing functions of the Brownian Motion at discrete times by 
functions of the diffusion at discrete times. Indeed, if we add the dw components 
Wt\ . . . , to the d components of the diffusion = {B]'\ . . . , 5f it IS easy 
to conclude that the diffusion obtained in this way generate the same filtration as 

0<t<T- 

Let f/' G L'^iJ^T] ^) and {X\ Z') the cJ- valued solution of {M + D) with = t/'. 
We are going to show that if ?7' — > f/ in L'^{J^t]'^), then (X') tends, for the distance 



to a process X ending at U and that there is a process Z such that {X, Z) solves 
the BSDE (M + D). 

Definition 4.1.1 Let T be the set of all terminal values U & uJ of processes X 
solutions of (M + D) (i.e. such that there is a process Z with {X,Z) solution 
of (M + D)) and S be the set of all these uJ-valued processes. According to the 
uniqueness part, to every U in T corresponds a unique process {Xt)t in S such that 
Xt = U . Hence we can define a mapping a with c{U) being this process : 



It is obvious that c is one-to-one and onto, and that is uniformly continuous 
for the above distances. Next we want to prove the uniform continuity of c. This is 
the aim of the following lemma, whose proof is a (slightly) modified version of the 
one of theorem (5.5) from 

Lemma 4.1.2 c is uniformly continuous for the distances (5'-^-' and 5^'^\ 

Proof. Suppose that we are given two terminal values Ui and corresponding to 
two solutions (X^, Zt) and (X^', Z[). Then by Holder's inequality, 




c: r ^ S 

U ^ iXt)t. 

We endow the space T with S^^^ and S with S^'^\ the distances just defined. 




(4.1) 



and since ^ 6^ 




CE sup e 
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The constant C above is allowed to vary from one inequality to another, but it 
depends only on T, U and \1/ (but not on the processes X and X'). The second 
inequality is Doob's one applied to the submartingale (exp(At)\I'(Xi, X/))^; the 
third one is Cauchy-Schwarz's one and the last one uses the classical inequality 
X < rjx^ + and the boundedness of 6 which implies <C5'^. 
Inequalities (j4.H) and (j4.2|) together give : 

<5(2)((Xi), (X;)) < CE (^e^.^+^i'o"(ll^^ll'+ll^^ll')"^)'^ (5«(?7i, U^)^^ . (4.3) 

Now we specialize to each case : if / does not depend on take /z = 0; then 
(e'^*\l/(Xt, X^'))^ is indeed a submartingale, and with 77 = = 0, the inequality 
()4.3|1 gives the uniform continuity required. 

In the other case, (exp(y4f)\I^(Xt, X/))j is yet a submartingale with = 
moreover, taking in (j4.3|) r] = a (with a as in Corollarv l3.4.4j) leads to the conclusion 
again. □ 

Using the completeness of the space of all aJ-valued processes endowed with the 
distance 5^'^\ we get the following result as an easy consequence of the preceding 
lemma. 

Proposition 4.1.3 Let {U^)i he a sequence in T converging to U & L'^{J-'t',uj) for 
the distance 6^^^ and X' = c(U^). Then there is a (continuous) process (Xi)ig[o;T] 

such that 5^'^\X\X) 0. In particular, we have Xt = U a.s. 

The result of this part is completed by proving the next proposition. 

Proposition 4.1.4 There is a process Z such that the pair (X, Z), with X defined 
in Proposition \4.1.3^ is a solution of (M + D) with Xt = U , limit of the sequence 

Proof. First, it is clear from the definition of S^"^^ that X is cJ- valued. We give 
then two proofs, one for the 2;-independent case and one for the general case in 
nonpositive curvatures. 

In the 2;-independent case, we use a proof very similar to the one of Theorem 
(4.43) of jTU] based on the characterization of Proposition 12.4.11 : first we localize 
the processes (in order to work in the Op defined in Proposition 12.4.11) and then it 
suffices to notice that the submartingale property passes through the limit in /. 

In the other case, it is not possible to apply the preceding proof, as the drift 
depends also on z. In fact, we are going to find explicitly a process Z such that 
(X, Z) is a solution of (M + D). We recall that for each /, we have in M" (remember 
that [ZlY is the J*'* row of the matrix Zl) 

dxl = z\dw, + (-|r,,(xl) ( [zl] ' I [z^^ + f{Bl XI ZD) dt 
x^ = uK 
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First step : Find a process Z, limit in the sense of the processes Z'-. 
Let 

X''™ = (X',X") and Z''" = 

apply Ito's formula ()2.H) to ^(X''™); then use ()3.10|) to bound below the Hessian 
term and (j3.1ip to bound above the term involving /. Taking the expectation gives 




-E 



ZT - z\ 



ds < CE (^sup5'(X^,X; 

+CE r S\Xi,X:^)i\\Z%. + \\Znr)ds. 
Jo 



We know that the first expectation on the right tends to zero as / and m tend to 
+00. Moreover, if Ji denotes the last integral, then using a a as in Corollarv 13.4.41 
and the fact that 6^ is bounded above by on the compact set cJ x cJ, we get 



Ji < V2E ( sup5^(X^,X^ 
"2 



E 



T 



7I \ l2 I II rymll 2 



zTVMs 



< ^/^EUup5^(X^,Xr) ) E{e^Ioi\\^i\\'r+\\zn'r)dsY 



c 



< J-E(sup5^(X^,Xr 

\ Oi 



(4.4) 



As C does not depend on /, m 

E 







f 


II ZT - z\ 


h 





as — > 0. 



Now we use inequality (j2.3p to have a bound on the Euclidean norm 

2 



(4.5) 



E 



zT-zifds < c 



E 



T 



ZT - z\ 



ds 



+ E r 5\XlXT){\\Z\\\l + \\Z, 
Jo 



mu2 
s II r 



)ds 



The first term on the right tends to by ()4.5p . and an argument similar to ()4.4|) 
would show that it also holds for the second term. Hence 



E 



7m '7l\\2 



zTds 



Z,m— ►oo 



0. 



Now by completeness of the space L^([0; T] x f2), we have the required result : 

3{Z,)teL\[0-T]xa):{Zl)^{Z,). 
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Second step : {X, Z) is indeed a solution of equation (M + D) with terminal 
value U. 

In view of this, let us show that the following expectation tends to zero as / tends 
to oo : 



E 



- + ZidWs + J^^ (^-lT,,{Xi){[Zi]' \[ZlY 



Obviously, this expectation is bounded above by 
E(|f/-f/f)^ + \\Zl- Z.fds^ ' 



+ \T,,iXs)-T,,iXl)\.\i[Z,f\[Z4^)\ds^ 
+ E \T,,{Xi)\ \{[Z,f I [Z.n - {[Zi]'\ [Ziy)\ds^ 



i-T 

+ E ' 



^ \fiBix,,z,)-f{By,xlzi)\dsy 



We know that the first two expectations tend to zero; the third term tends to zero 
by dominated convergence (at least for a subsequence of (X'), but it doesn't matter 
since (X') is Cauchy). Let Ei denote the next expectation; then we can write 

El < ceQ^^ \\Zi-ZM\Zi\\ + \\Zs\\)ds^ 

< v^CE \\Zl - Zsfds^ ' E (^j\\\Zif + \\Zsf)ds^ ' • 

The first integral tends to zero and the second is bounded because Z'' converges in 
L^; hence Ei tends to zero. 

Finally, let E2 denote the last integral and use inequality ()2.5p to obtain 

< E {S{XlXs){l + \\Zl\\ + \\Z,\\) + \\Zl -Z,\\) ds^ ; 

using an argument similar to ()4.4|) . this quantity tends also to zero as / tends to 00. 
Hence the limit in of Xl is 



U-j^ ZdWs- (^-]^T,u{Xs){[Z,t\[Zs]') + f{Bl,Xs,Zs)^ds- 



We know that it is also Xj, so by continuity 

-T r-T 



a.s.,\ft, Xt = U- ZsdWs-J^ (^-hjj,iX,)i[Z,f\[Zsy) + f{By,X,,Zs)^ds. 
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That finishes the proof of the proposition for a general / in nonpositive curvatures. 
□ 

As a consequence of Propositions 14.1.31 and 14.1.41 it suffices to work with a 
dense subset of L'^{J-'t',uj), i.e. to consider terminal values U that are written U = 
G{Bl,Bl, ...,Bl) where G G C^{W^,uJ) (in fact, we can obviously take tg = T). 

A further step of simplification is possible (see [IS]): conditioning by J-'t^^i and 
working over the time interval [tg_i;T], B^^, . . . ,B^^_^ can be treated as constant. 
Then, if we know a solution (X^, Zt) to (M + D) on [tq_i;T] with Xt = G{BJj.), 
the problem is to reach But this variable is in fact a measurable function of 

Sf^, Bf^, . . . , Bf^_^. A density argument enables us to suppose this function smooth 
and constant off a compact set. Thus by induction the problem is solved if we can 
find a solution to (M + D) with terminal value F{Bl) with F e C^{W^,uj). 

So to prove the existence result for any oJ-valued terminal variable it suffices 
to solve equation [M + D) with a terminal value U that can be written F{B^) where 
F e C^(]R'^,cJ). This is the aim of the next paragraphs. 



4.2 Approximation by BSDEs with Lipschitz coefficients 

We deal with a BSDE whose coefficients Tjklx) (we recall that it is a vector in R"; 
see the introduction) and f{b,x) (or f{b,x,z)) are defined only for x in the open 
subset O of M", and with a quadratic term in Zt; but we would like to apply the 
existence result of [21] to BSDEs with Lipschitz coefficients. Our purpose in this 
part is to define such a BSDE, derived from the initial one and defined on all M". 

Firstly, we extend the definition of the BSDE (M + D) to the whole space M" : 
let be a smooth function on M" with compact support in O and such that 0=1 
on uj. We will explicit a little more later in Subsection 14.41 Then we extend 
the drift / and Christoffel symbols to the whole space M" by letting , for b G M*^, 
fib,x) = (j){x)f{b,x) (or f{b,x,z) = (l){x)f{b,x,z)) and f{x) = (p{x)T{x). 
The new BSDE defined on all (we keep X and Z for the notations) is : 

it,,{Xt){[Ztf I [Ztf) + f{BlXt, Zt)) dt 




Remark : We write f{Bt,Xt,Zt) but when it is not pointed out, it should also 
be interpreted as f{Bt,Xt) as well. 

Now let e €]0; 1[ and : R+ — > R+ be a smooth nondecreasing function such 
that Si;{t) = iff t G [0; ^], and is linear for t large enough. Then, for z G R"'^^, 
we define hgi^z) = S£{\\z\\) and z = Yq^Xz)' have the following 

Lemma 4.2.1 The functions /i^ and z defined above satisfy the next assertions. 

(i) If \\z\\ < ^, then h^^z) = and z = z; 

(a) z ^ z is smooth, bounded and Lipschitz on R"*^^; 

(Hi) {x,z) g{x,z) := Tjk{x){[z]^\[zy) is smooth and Lipschitz on R" x R"'^'*'. 
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(iv) {b,x,z) I— > f{b,x,z) is bounded and Lipschitz on R'^ x M" x M"*^^. 
As a consequence, the function 



7 : 



^'^^ X M" X M""*^ 



{b,x,z) 



h,,{x){[zf\[zy)-fib,x,z) 



is bounded and Lipschitz. 



Proof. 

(i) It is obvious on the definition. 

(ii) The smoothness and boundedness of z are clear. Now we just prove that 
11^ ~ z'W < C\\z — z'W, considering the 3 cases below. 

• \\z\\ < ^ and \\z'\\ < ^ 
Then \\z — z'\\ = \\z — z'\\. 

• \\z\\ > - and \\z'\\ > - 

Then, as 1 is bounded and /i^ Lipschitz, we can write 



z — z' 




1 + h,{z) 


+ \w\\ 



\z-z'\\ < 

< ||z - ^'11 + 

< C\\z-z'\\. 



1 



1 



l + h,iz) l + h{z') 
K{z') - h,{z) 



l + hJz) 



• ||z|| > ^ and \\z'\\ < ^ 
Let z" G [z; z'\ with ||2;"|| = -\ then, using the first two cases, we have 



\z — z 



< 



\z~z"\\ + |U' < C{\\z-z"\\ + -^"11) = C\\z-z\\. 



(Hi) We have 



\gix,z)-gix\z')\ < |f,,(x) - f,fc(x')| ■ 

+ |f,,(x')|-|([^fim-([^1[^^')l- 

There is no problem with the first term on the right because Tjk is smooth on 
with compact support and z is bounded. Using (ii), it is easy to deduce that the 
second term is bounded above by C\\z — z'\\ and that finishes the proof. 
[iv) Recall that for x,x' G O, 

\fib,x,z)- fib',x',7)\ < L'{{\b-b'\ + + ll^ll + ||7||) + \\z-^\\) . 

If Os denotes the compact support of 0, this inequality and imply that f{b, x, z) 
is bounded, uniformly in x G O^, 6 G M*^ and z G W^'^^ . Thus / is bounded on 
M*^ X X because for x outside Os, / = 0. Besides, we consider the 3 

following cases to study the Lipschitz property : 
• x^x' ^ Os 

Then f{b,x,z) — f{b',x',z') = and it is obvious. 
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• X e Os,x' ^ Os 
Then, as is Lipschitz, 



\f{b,x,z)-f{b',x',z')\ = \(j){x)f{b,x,z)\ 

= \^{x)-4>{x')\-\f{b,x,z)\ 
< a\x — x'\. 

Then 

\fib,x,z)-f{b',x',^)\ = \(f){x)f{b,x,z)-(f){x')f{b',x',7)\ 

< \4>{x)-^{x')\-\f{b,x,z)\ 

+ \^ix')\-\fib,x,z)-f{b',x',7)\ 

< a\x — x'\ + f3i\\'z — z'W + '~f\b — b'\ 

< a\x — x'\ + l32\\z — z'W + '~f\b — b'\ 

since z \-^z is Lipschitz by (ii). The proof of (iv) is completed. □ 
Now we can introduce a new BSDE 

The interest of this new equation hes in the following result. 



dX^ = ZtdWt-j{BlXf,Zf)dt 

= u. 



Proposition 4.2.2 The equation (M + D)^ has Lipschitz coefficients; it has a 
unique solution (Xf , Zf) e M" x M"'^^- ^^^./^ ^/^^^ 

E f sup \X^\A < oo and [ \\Z^fdt^ < oo. 

\te[0;T] J \Jo J 

Proof. The first part is a consequence of Lemma [4.2. 11 then existence and unique- 
ness are classical results of □ 



In Subsections 14 . H| lOl and below. / is also supposed to be a function such 
that {b,x,z) 1—^ f{b,x,l) and 7 are functions with all their partial derivatives of 
order 1, 2 and 3 bounded. 

4.3 Existence of a solution of (M + D) on a small time in- 
terval 

In this paragraph we consider a terminal value U = F{Blj^) {F G C^{M.'^,uj)) which 

lies in uJ and we show that for an e small enough, the solution {X^, Z[) of (M + D)^ 

is also a solution of the BSDE (M + D) on a small time interval [Tf ; T]. This result 
relies on Proposition 14.3.21 which gives the very strong condition that Z is bounded 
a.s. 

Firstly, let us give a classical link between BSDEs and PDEs. 
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Proposition 4.3.1 Let e g]0; 1[; we use the notations introduced in paragraph \4 . ij| 
(recall in particular that the function 7 depends on e). 

Consider for Us = {ul, . . . , m") : [0; T] x M'^ — > M" the following system of quasilinear 
parabolic partial differential equations 

= Cus{t,x) +'y{x,u^{t,x),{\/a,Uea){t,x)) ^^^^ 



UeiO,x) = F{x) 

where VxUe is the n x d matrix whose rows are {'VxUl)i=i^,,,^n, the partial derivatives 
of the components of Ue with respect to space; moreover, 

/:u, = {Lul,...,Lu^) (4.7) 

is a vector in M", and 

t,j=l 1=1 

is the infinitesimal generator of the diffusion ()1.2j) . 
Then 

(i) This equation has a unique solution Ue in C"'^'^([0;T] x M'^jIR") (i.e. has 
continuous first derivative with respect to time and continuous second derivatives 
with respect to space); 

(ii) A.S., XI = u,{T- t, Bf) Vt; 

(ill) Vt, a.s., Zl = V^UeiT - t, By)a{By). 



Proof. The assertion (z) is a classical result; for a probabilistic proof, see Theorem 
3.2 of [SHj, and for a proof in dimension 1 (i.e. n = 1) with less regular functions 7, 
see (21. 

Then (ii) and {Hi) follow easily : both Ito's formula and equation ()4.fjj) lead to 

dueiT - t, Bl) = -7(5f , Ue, iV.,Uea)iB!))dt + {\/.Uea)dWt 
uMB't) = F{Bl) = U. 

Thus 

(n,(T - t, Bt), V^UeiT - t, Bt)a{By)\ 

is a solution of BSDE (M + D)e. But (X*^, Z'^) is also a solution and this BSDE has 
Lipschitz coefficients, so it has a unique solution. Then {ii) and {Hi) follow. □ 

Remark. Since V is continuous, the trajectories t ^— > Zl can be taken 
continuous. 

Now we give the main result of this paragraph. 
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Proposition 4.3.2 There is an e g]0; 1[ and a Tf G [0;T[ (depending on e) such 
that a.s. for t G [Tf;T], we have \\Zf\\ < K This means, with the notations of 
paragraph \4-^ that hs^Z^) = so = Zf and {Xf^Zf) is a solution of BSDE 
(M + D) on the time interval [Tf; T]. 

Proof. We only deal with a drift / depending both on x and z, the 2;-independent 
case being similar (and easier). 

We let be the global Lipschitz constant of the function 7 (independent of h): 

V(x,;z),(a;',z') G M" x M"^^, 

\l{h,x,z) - 7(6,0;', < ae{\x-x'\ + ll^-^ll). (4.9) 

(in general, will tend to 00 as e decreases to zero). 

We have to show that \\Zl\\ = WVxUeiT — t, i?f)cr(_Bf )|| < ^ a.s., so it suffices to 
check that Ue is -Lipschitz with respect to the space variable (||o"||^ denoting 

the supremum of ||o"(6)|| for b G W^). This work will be achieved in several steps. 

First Step : We consider two solutions (Y, Z) and (Y, Z) of (M + D)^ corre- 
sponding to two terminal values Yt and Yt. We set 6Yt = Yt~Yt and 6Zt = Zf — Zt. 
Then 

^Yt + j^ 5ZsdWs = 5YT + (^jiBy,Y„Zs)-^iBy,t,Zs))ds. 
Independence gives : 



T 



El^FtP^E / \\6Zsfds = E 



T 



6Yt+ {^{By,Y,,Z,)--f{By,Y,,Zs))ds 



2^ 



Using Holder's inequality, we obtain : 

E\6Yt\^ + e(^J^ \\6Zsfds'^ < 2 E\6Yt\ 



+ 2(T-t)E| // \j{By,Y„Z,) 
-j{By,Ys,Zs)?ds 



Hence using (j4.9p we have 



+ 4a2(T-t)^ E(/f 

+E (^J^^ \\6Zsfds 
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Second Step : We want to remove the terms E ^ to get an equa- 

tion with Y only; this leads to make the following assumption (which will be true 
in the end of the proof) : t G [Tf ; T] with Tf G [0; T[ such that Aal{T - Tf ) < 1. 

Then, if t G [Tf ; T], we get 

E|(5yt|2 < 2E|(5FTr + E \5Ys\^ds 
And Gronwall's lemma gives 

'E\5Yt\^ < 2e^-*E|5yTp- (4.10) 

Third Step : Let us choose Yt = F{B^/) and Yt = F{B!f) (where (5*'^) 
denotes the diffusion starting at x at time t); then 

5Yt = u,{T -t,x)-u,{T -t,x) (4.11) 

and 

El^lrl' < T|E - E^'p) < - x\\ (4.12) 

with Li? the Lipschitz constant of F, and the last inequality is a well-known result 
of L^-continuity with respect to initial conditions. Hence from ()4.10|) 



|Me(T — t, x) — Me(T — t,x)\ < Lf\ ^Ca^hG^ ^\x — X 



and 



\V^Ue{T -t,x)\\ < L^J2a,f,e^ 



(In fact, the right hand side is given up to a positive multiplicative constant, due to 
the different norms used; it doesn't matter in the sequel). 

Conclusion : Let e g]0; 1[ such that ^ > ^a^^oL f \ f^^C^^ifi^ . Then the preceding 
proof shows that if Tf G [0; T[ is such that Aol1{T — Tf ) < 1, then for t G [Tf ; T] we 
have 

a.s. ||Zf|| = ||V,.n,(T-t,i?f)a(5f)||<l. (4.13) 

In fact, using the continuity in t of Zf, this inequality holds for any t, a.s.; that 
finishes the proof of the proposition. □ 



4.4 A solution of (M + D) on a small time interval 

The framework is the same as the one introduced in Subsections 14.21 and 14.31 The 
aim of this section is to prove that the above solution (Xj, Zt)t (we omit in this 

paragraph the superscript e for notational convenience) of (M -|- D) is, under an 
additional condition on the drift /, a solution of (M + D) (in fact we prove that 
(Xt) remains in aJ, the compact where the terminal value U lies). We recall that 
ZU = {x < c} is the sublevel set of a smooth convex (for the connection F) function x. 
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defined on an open set O relatively compact in M", and moreover that uj is relatively 
compact in O. We make the following hypothesis 

{Hg) f is pointing strictly outward on the boundary dZJ of to. 
It means that 

V(6, x,z) : xe duo, inf {Dx{x)\f{h, x, z))., > C > 0, (4.14) 

b,x,z 

where (■|-)r denotes the Riemannian metric tensor (if /-*- is the component of / 
orthogonal to duo = {x = c}, it is equivalent to require that inf {Dx{x)\f-^{h, x, z))r 

b,x,z 

be bounded below by a positive constant, since x is constant (equal to c) on duJ). 

Remark : This condition arises naturally in the deterministic version of equa- 
tion (M + D) 

( dxt = f{xt)dt 
\ Xt = u 

where u is deterministic and so Zt = for any t (in fact, as we shall see in Sub- 
section 14.61 the natural condition is to require that the infimum in (|4.14p be only 
nonnegative). 

Proposition 4.4.1 Suppose that x is strictly convex on O (this means that Hess x 
is positive definite). Then, under the assumption (Hs), the process {Xt)t£[T~-T] 
mains in uJ, i.e. {Xf, Zt) is a solution of (M -|- D) on the time interval [T-f; T]. 

Proof. Once again, the goal is to construct a nonnegative submartingale, null at 
time T, which vanishes if and only if the process X is in uJ. The proof will be split 
into three steps. 

First Step : Framework of the proof. 

Suppose that c > and x reaches its infimum at p E u with xip) = 0. Then con- 
sider the following mapping, defined on a normal open neighbourhood of uJ centered 
at p (so it is in particular a neighbourhood of in M") 

y ^ I =i 

VC2||yP-Fc-x(y) 

for C2 > small enough, it is a diffeomorphism from an open set 0\ (relatively 
compact in O and containing uo) onto an open neighbourhood of 5(0, 1), such that 
To is sent onto 5(0, 1) (in fact, it is sufficient to take C2 < \^xi "where A^^ denotes the 
(positive) infimum on 0\ of the eingenvalues of Hess x)- Using this diffeomorphism, 
we can work in a local chart O (take O := A^) such that cJ = -8(0, 1) and x(0) = 0. 

Second Step : Construction of a "nearly convex" function H. 

Choose p > 1 such that B{0,p) C O; then we take for the mapping (/) (see the 
beginning of Subsection 14. 2|) a smooth function M" [0; 1] equal to 1 on B{0, 1) and 
to outside B{0,p), which has moreover spherical symmetry. This gives T = (pT 
and f = (f)f. Then define on 

k{x) = 0(x)x(x) -l- (1 — (f>{x))a{x) 
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where a{x) = a\x\ and a > is chosen so that k{x) < c if and only if x G 5(0, 1) 
(take for instance a > sup^j-g.p) x)j k is clearly a nonnegative smooth mapping. We 
would like to have a mapping which vanishes on uJ, so we let H = h o k where 
h : ]R+ ]R+ is a smooth convex (so nondecreasing) function, vanishing on the 
interval [0; c] (only) and growing linearly at infinity. 

The mapping H so defined is convex for the connection F on B{0, 1) and outside 
B{0,p). Indeed, on B{0, 1), H = and outside B{0,p), H = ho a which is convex 
for the flat connection (= F). 

Third Step : We work on the time interval [T^^;T]. The aim of this step is 
to show that, for A large enough, the process {e^^H{Xt))t is a real submartingale. 
Reasoning as in the uniqueness part, we apply Ito's formula : 

e^*ff(XO - e^^f i/(XTf ) = / e^'DH{X,){ZsdW,) 

[j2'['ZsY^ HiXs)['Z,Yj ds 
+ / e'WH{Xs) ■ f{By, X„ Zs)ds 
+ [ Xe^'H{Xs)ds. 

The stochastic integral is a martingale because DH is bounded; it remains to 
prove that the bounded variation term is an increasing process, i.e. show the non- 
negativity of the sum 

-J2'['^y'^H{x)['zf + DH{x)f{b,x,z) + XH{x). 

i=l 

But Hess H{x) > h' {k{x))B.ess k{x) since k is nondecreasing (see for instance 
(4.36) p 42 in [ID]); then it suffices to prove the nonnegativity of 



(1 \ 
-J2'['^f^ + Dk{x)f{b,x,z)] +XH{x). (4.15) 

Remark that h'{k{x)) > and that, because of the boundedness of the process 
Zt (according to ()4.13j) ). it is sufficient to consider z such that ||z|| < -. 

Let Ai, A2 and A3 denote the three terms in this order in the sum ()4.15|) . 

• if a; G -8(0, 1), Ai and A2 vanish because if = 0, so the sum is nonnegative; 

• if x G ^-8(0, p), Ai > because k = a is convex for the fiat connection and ^2 = 
because / = 0; then the sum is nonnegative. 

• if X G -8(0, p) \ -8(0, 1), we consider two situations : firstly when x belongs to a 
neighbourhood of the sphere S{0, 1) (i.e. when \x\ g]1; 1 + ri[ where > is to be 
determined). Using continuity, Hess k{x) is nonnegative (in the sense of matrices) on 
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]1; 1 + (with rji > sufficiently small) because if |x| = 1, Hess k{x) = Hess x{^) 
and X is strictly convex on -8(0, 1) for the connection F. This gives Ai > 0. 
For A2, using hypothesis (Hg), we get that for x G duJ and ||2;|| < i, 

Dk{x) ■ f{b, X, z) = Dx{x) ■ f{b, x, z) = Dx{x) ■ f{b, x,z)>C>0. 

So by uniform continuity and hypothesis ()1.4|1 and ()1.5|) . there is an 772 > such 
that 

Wb, \/y:l<\y\<l+V2, : ||^|| < ^, Dkiy) ■ fib,y, z) > 0. 

In particular, A2 > for 1 < |x| < 1 + 772- Let rj = min{rii,T]2, p — 1). Then 
+ ^2 > if 1 < |x| < 1 + < p, and the sum (j4.15|) is nonnegative. 
On the other hand, if |x| e]1 + rj; p[ (i.e. if |x| is "far" from 1), then Ai + A2 is 

bounded above since H is smooth, and f{b,x, z) is bounded since z is bounded. But 

we have constructed H so that H{x) > > if |a;| > 1 + r^. Then we can choose 

A > such that Ai + A2 + A^ is nonnegative. 

Conclusion : The end of the proof now goes on by classical arguments : the 

process {e^^H{Xt))t£[Tf;T] is a nonnegative submartingale for the A chosen above, 

null at time T. Then Vt G [Tf;T], H{Xt) = 0; and the mapping H has precisely 

been chosen so that this implies G oJ. □ 

Remark : In our local coordinates, let Drx{x) denote the radial derivative of x 
at X and the radial component of /. Then, since x = c on the sphere S{0, 1) = dU, 
we have 

Vx G duJ, Dx{x) ■ fib, x, z) = Drxix) ■ f{b, x, z). 

But if Hess x is positive definite, Drx{x) is a positive real number and in our coor- 
dinates, {Hg) is equivalent to require that /'"(&, x,z) be bounded below by a positive 
constant, independent of x G du and of b, z. 

Now we come back to the end of the proof of the existence. 

4.5 The solution on the whole interval [0; T] 

According to Subsections 14.31 and 14. 4t (Xf, Z^) is a solution of BSDE (M + D) on 
the time interval [Tl;T]] moreover, e and must verify (see the conclusion in the 
proof of Proposition I4.3.2p 

f 7 > |k||oo^Fv/2CV^ 

\ 4a2(T-Tf) < 1 

where tends a priori to 00 when e goes to zero. The aim of this subsection is to 
show that it is a solution on the time interval [0; T]. 

In general, if we take back the proof of Proposition I4.3.2| with now Tf and 
Us{T — Tf , Bl^e) as terminal time and variable, we get a solution of (M + D) on a 
time interval [T|; Tf]; but it is easy to show that e and T2 must verify now 

f i > (||a|UL^v/2a;^)v^ 

1 4a2(Tf-T|) < 1; 
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as a consequence, the length of the interval [T|;Tf] may be less than the one of 
[Tf ; T] and repeating this method inductively could lead to a solution on an interval 
]To;T[ with Tq > only (this means that the solution explodes). 

In fact, the existence of \I' prevents the solution from exploding and allows to 
build inductively a solution of (M + D) on [Tf ; T], [T - 2{T - T^); T], [T - 3(T - 
Tl);T], ... and so to get at the end a solution on [0;T]. 

Proposition 4.5.1 Suppose that assumption (Hs) holds and that x is strictly con- 
vex on O. If e > is small enough, then {X^, Z^), the solution of BSDE (M + D)^, 
is also a solution of BSDE (M + D) on the whole time interval [0;T]. 

Proof. Consider Tf and Me(T — Tf,i?^e:) as terminal time and variable; applying 
the two preceding sections, we get a solution of (M + D) (with Xf e aJ) on a time 
interval [T|; Tf] where e and T| verify 

J ~ ^ ||o"||oo-^U£ 

\ 4a2(Tf-T|) < 1 

{Lu^ is the Lipschitz constant of Ug for the space variable, uniformly in t on [T|; T]). 
As in Subsection 14.31 for t G [T|; T[ we let {Yg, Zs)se[t;T] and (1^, Zs)s<^it-T] be two 

solutions of BSDE {aF+D), such that Yt = F{B^/) and Yt = F{B^/) (by which 
we denote diffusions starting at x or f at time t). According to Subsection 14. 4| these 
two processes remain in ZJ, so we can make use of the function Then the same 
inequahties as in ()4.2|) give 




(the second inequality is obtained by letting r/ = in the ^-independent case; in the 
other case, by letting r/ = a as in Corollarv 13.4.41 and p = 2). 

The equivalence of the Riemannian and Euclidean distances on EJxcJ and \& ~ S^, 
together with pTT|l and give 

lusiT — t,x) — Ue(T — t,x)\ < CoLp\x — x\, 

where Cq depends only on the diffusion {B^)t, T, U and the drift / (then L^^ = 
LpCo). 

Consequently, if we take e such that 

- > Iklloo^FCo, (4.16) 

e 

Tf s.t. 4a^(T — Tf) < 1 and T2 s.t. ia'^iTl — T2) < 1, we get a solution on an 
interval [T| — Tf] with the same length as [Tf;T]. Repeating the same method 
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inductively with the same e at each step, we obtain that {Xf,Zf) is a solution of 
(M + D) on [Tf ; T], [T - 2(T - Tf ); T], [T - 3(T - Tf ); T], ... and so on the whole 
interval [0;T]. □ 

Remark : As = VxUe{T — t, B^)a{B^), it is a straightforward consequence 
of the proof that a.s., for all t, \\Zf\\ is bounded above by -. 
If we sum up the results obtained, we have the 

Proposition 4.5.2 We consider BSDE (M + D) with a terminal value UinuJ = 
{x ^ c}. Suppose that f is a function which verifies conditions p.4|) . (jl.Sp and 
{Hs). If moreover X is strictly convex (i.e. Hess x is positive definite), then 

(i) If f does not depend on z, the BSDE has a (unique) solution such that X remains 
in uj. 

(ii) If M is a Cartan-Hadamard manifold and the Levi-Civita connection is used, 
then the BSDE has a (unique) solution such that X remains in ZJ too. 

The last paragraph is devoted to generalize Proposition I4.5.1?! to drifts / which 
are less regular and satisfy a weaker hypothesis than {Hg). 

4.6 The solution for general / 

Let / be a function verifying ()1.4|1 . and the hypothesis 

(H) f is pointing outward on the boundary of uJ 

introduced in Subsection 11.41 This means that 

y{b,x,z) : X e diu, {Dx{x)\f{b, x, z))r > 0; 

equivalently, if /-*- is the component of / orthogonal to dZJ = {x = c}, we may 
require {Dx{x)\f^{b,x, z))r to be nonnegative. This condition is obviously weaker 
than {Hg). 

Firstly, we will derive from equation (M + D) new BSDEs, each one having a 
unique tJ-valued solution. Then we will show that these solutions converge to the 
solution of (M + D). The function x will supposed to be strictly convex and calculus 
will be achieved for a drift / depending both on x and z (the case when / does not 
depend on z being simpler). 

Let us consider the local chart O defined in the first step of the proof of Propo- 
sition l^7¥?Tl Remember that aJ = -6(0, 1) in these local coordinates and remark that 
hypothesis (H) means that the radial component f^{b,x,z) of f{b,x,z) is nonnega- 
tive for X G duJ (see the remark at the end of Subsection I4.4|l . Extend the mapping / 
to M"^ X M" X W^"^ by putting f{b,x,z)=0ifx^O and define (on M'^ x x M'^'^^) 
the convolution product for / G N* = / * pi where pi{b, x, z) = lp{l\\ (6, x, z) ||) and 
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p : M+ ]R+ is a bump function (i.e. a smooth function with p'{0) = 0, p = 
outside [0; 1] and f^_^ p{u)du = 1). Besides, let us define a function gi by 

gi{b,x,z) = fi{b,x,z) + jx, (4.17) 

where A is a positive constant which will be chosen below; we introduce on O the 
following BSDE for F e C^iR"^) and ?7 G cJ 



(M 



/ dXt = ZtdWt+{-'jT,k{Xt){[Zt]''\[Zty)+giiBlXt,Zt))dt 
\Xt = F{B<.) = U. 



Lemma 4.6.1 Let Oi be an open set such that uJ G Oi G Oi G O . Then for I > Iq 
the functions gi are smooth on M*^ x Oi x M^'^w-' and verify a Lipschitz condition like 
(12. 5p with the same Lipschitz constant L" . Moreover, they also verify a boundedness 
condition like (|1.5p . 

Proof. For a; G Oi, we write 

fi{b, x,z) = / /((6, X, z) - {(3, y, w))pi{(3, y, w)d{p, y, w). 



Thus, as soon as dist{Oi,0) > I/Iq > l/l, pi{P,y,w) = if \y\ > dist{Oi,0); so 
the integrand vanishes if x — y ^ O and we can use the Lipschitz property (j2.5|) 
of / on O. Then the properties of convolution give the result for As it also 
holds obviously for the functions x ^ A/l x, we have the result for gi. The second 
assertion is an easy consequence of conditions ()1.4|) and (jl.Sp for the drift /. □ 

In the sequel, we will consider the sequence {gi)i only for I > Iq. 

Proposition 4.6.2 For every I, the BSDE {M + D)i has a (unique) W-valued solu- 
tion {X\ Z^); moreover, there is an e > 0, independent of I, such that a.s., \\Z[\\ < ^ 
for any t. 

Proof. Using Lemma I4.6.H we apply Subsections 14.21 and 14.31 to {M + D)i. We 
get an ei and a Tf' for every /; but the Lipschitz constant of 7 is independent of I 
(since the gi have the same one), so the proof of Proposition 14.3.^ shows that we 
can choose Si = e and Tf' = Tf , independently of /. 

In order to apply Subsection 14.41 we need to prove that gi verifies condition (i/^); 
in fact, we have seen (see in particular the remark at the end of Subsection 14. 4|) that 
it suffices to show that 

V(6,x,;z) : xGdu, Ikll < mf c/[(6, x, ;z) > C > 0, (4.18) 

£ b,x,z 

It is easy to see with the properties of convolution that, for a constant C depending 
only on e and /, 

V/, V6, Vx edLJ,^z: \\z\\ < -, \fi{b,x,z) - f{b,x,z)\ < y. 
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This and the nonnegativity of f^{b, x, z) for x G duo give 

V/, V6, Vx G du, Wz : < ^, f[ib,x,z) > -j. 

Now if we take in pTTjl A = C + 1, obviously holds with C = and 

Subsection 14.41 can be applied. 

Then the results of Subsection 14. 51 hold (in particular Prop osit ion 14 . 5 . 2]) with the 
same e for every I (this comes from ()4.16|) . remarking that the constant Co in this 
inequality is independent of /). The proof is completed. □ 

Now we prove that, as expected, the limit of these solutions when / —>■ oo solves 
equation (M + D). 

Proposition 4.6.3 For n G N and s G [0;T], we note for simplicity Vp : = 
{By,Xg,Z^), When I tends to infinity, the above solution {X\ Z'') converges (for 
the usual norms for X and Z) to a pair {X, Z) which solves equation (M + D). 

Proof. For / and m in N, we note 6Xt = X^ - , 6Zt = Z\ - and 

Ai^^ = e£ \SX,\ (^IfmiV^) - fiiV^)\ i-Xr ) ds. 

Applying Ito's formula (in M") to I^X^p between t and T, we get (note that a 
dot stands for the inner product in M") 



\SXt\^ = I 25X,-i5ZJWs) + 2SX, ■ (giiV^) - g^iVD) ds 

n 

26X, . {r,,{x%zmziy) - v ,^[x-)(\z^f\\z^\^)) ds 

\6Zsfds. (4.19) 

Using the uniform boundedness of the (Z') (proved in the preceding lemma), we 
can bound above the integral involving the Christoffel symbols by C \6Xs\'^ds; 
besides, for the second term on the right, we write 

5X, ■ (g^iV^) - gUVn)ds 



< Am + j^ \6Xs\-\fUV^)-fm{Vsn\ds 

< Ci(^Ai,m + j^ \6Xs\i\5X,\ + \SZs\)ds^ 



< C(Ai^„,+ / \6Xs\^ds) + 



\6ZJ\ ds 
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where C is independent of / and m. Then, we obtain by taking the expectation in 



^5Xt\'^ + \\6Z,fds<c(^J^ E\6Xs\'^ds + Ai^rn^ . (4.20) 

Gronwall's lemma gives E|5Xtp < CAi^^, where C is again independent of I and 
m. Moreover, using ()1.4|) . fi converges uniformly to / on M'' x cJ x -8(0, r) for any 
r > (with 5(0, r) = {z e M"''^ : < r}) and X\ X"" are bounded, so Ai^^ 
tends to zero when /, m tend to infinity; therefore {X'')i converges to a process X in 
L2(n X [0;T]). 

Using ()4.2()|1 again, we get 

E / \\5Zsfds < C\^- 
Jo 

hence the sequence of processes (Z') has also a limit in L'^{Q x [0; T]); let Z denote 
this limit process. 

The pair (X, Z) solves BSDE (M + D) and X is tJ-valued; the proof is just 
an adaptation of the second step in the proof of Proposition 14.1.41 This remark 
completes the proof. □ 

Remark : As a consequence, a.s., \\Zt\\ < ^ for any t. 

According to Subsection 14.11 this result can be extended to every cJ- valued and 
jFT^-measurable terminal variable U. Then Theorem 11.4.11 of existence and unique- 
ness of a solution follows. 

Note that uniqueness and existence hold in particular on any regular geodesic 
ball (or geodesic ball if the sectional curvatures are nonpositive) . 



5 Applications and related PDEs 
5.1 The martingale case 

The drift / = verifies hypothesis (H). Hence in this case the results of this paper 
apply to the martingale case. As already underlined, any regular geodesic ball veri- 
fies the condition of Theorem ll.4.11 so we recover the well-known results of existence 
and uniqueness of a martingale with prescribed terminal value in such domains (see 
[T3]). These results hold in nonpositive curvatures, they will be achieved in positive 
curvatures elsewhere. 



5.2 The one-dimensional case 

The nonpositive curvature case gives the existence and uniqueness of a solution to 
the one-dimensional BSDE 



(Eh 



dXt = ZtdWt - T{Xt)Z^ + f{Bl Xt, Zt)dt 
Xt = U 
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for a bounded terminal condition U, a drift / satisfying (jl.Sp . ()2.5|) and any smooth 
function F defined on M. 

Note that a change of coordinates (in fact a reparametrization of the one- 
dimensional manifold by arclength) reduces equation {E)i to 

r dXt = ZtdWt + f{BlXt,Zt)dt 
\Xt = U; 

moreover, it is a very particular case of the results of Kobylansky in . 

One can ask whether such results can be extended to higher dimensions. In fact, 
the original problem is geometric and to deal with general BSDEs, we would start 
with smooth functions (^ij) and should give conditions in order to interpret these 
functions as the Christoffel symbols of a given Levi-Civita connection. This problem 
is out of the scope of this paper. 

5.3 Case of a random terminal time 

In this paragraph, we will only sketch the proofs. 
We are interested in the following equation 

dXt = ZtdWt + {-\T^k{Xt){[znztY) + f{Bl Xt, Zt)) dt 

Xr = 

where r is a stopping time with respect to the filtration used and f/^ is a tJ-valued, 
Jv-measurable random variable. It is the counterpart of equation (M + D) on the 
random interval [0; r]. 

Let us first consider the case of a bounded stopping time r, i.e. t < T where T 
is a deterministic constant. We have the following result : 

Theorem 5.3.1 We consider BSDE (M + D)^. with To = {x < c\ and t <T a.s. 
If f verifies conditions ()1.4|) . ()1.5p and (H), and if x is strictly convex (i.e. Hess x 
is positive definite), then this BSDE has a unique solution {X,Z), with X G tJ, m 
the same two cases as in Theorem \1.4.1\ 

Proof. First remark that the uniqueness part goes the same as in the deterministic 
case; for the existence part, it can be completed in several steps : 

First Step : We work in local coordinates introduced in Subsection 14.41 Let 
ci > c and put tJi = {x < ci}; suppose that ci is such that tJ C cJi C O and 
that X is yet strictly convex on uji. Now let be a cut-off function with 0=1 
on uJ and = outside Wi. For any nonzero integer /, we define a new drift 
by fi{b,x,z) := (j){x){f{b,x,z) + (l//)x); note that /; verifies hypothesis (H) with 
respect to aJi and {Hg) with respect to aJ. 

Then solve path by path on [r; T] the following differential equation 

r dXl = fi{BlXl,0)dt 

\Xl = U^ l^.lj 



{M + D)r 
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and set = + J' fi{B^ , Xl,0)dt. Since /; vanishes outside lui, U is in ui; 
besides, f/' is jF-r-measurable. 

Second Step : Considering the random variable f/' and the drift fi introduced 
in the first step, we solve on [0; T] equation (M + D) with drift // and terminal 
value Xlp = f/'. The hypothesis of Theorem 11.4.11 are satisfied considering the set 
ZJi instead of ZJ. So this BSDE has a solution (X', Z') with X' G tJi. 

We condition by !Fr and consider the above solution (X', Z'-) on the random time 
interval [r; T]. It is a solution on this interval of BSDE (M + D) with drift // and 
terminal value f/'. The uniqueness property for such equations implies that is 
the solution of equation (jSHJ, i-e. 

Vt G [r; T], Xl = U' - fiiBl X^, 0)ds 
and Zl = for T < t < T. In particular, 

xl = u'- £ fi{By,xlo)ds = u\ 

We now show that actually, (X')o<t<,- remains in uJ, and not only in tJi. For this 
purpose, we want to construct as in Subsection 14.41 a submartingale which is written 
as (e'''*if(Xj')), where now H = ho x, since X' is cJi-valued. Recall from ()4.15j) that 
the keypoint is to show for x G cJi the nonnegativity of 

d-w \ 

-J2'['zYRess x{x)['zY + Dx{x)Mb,x,z)j + XH{x). 

In fact, as x is strictly convex on the compact domain cJi, we have Hess x ^ otid 
(in the sense of matrices) for a > and it turns out that it suffices to prove the 
nonnegativity of 

h'ixix)) Qalkf + Dxix)fi{b, X, z)^ + XH{x). (5.2) 
From p.4|) and ()1.5p we deduce 

\Dx{x)fi{b,x,z)\ < Cil + \\z\\) <^a\\zf, 

the last inequality holding for \\z\\ large enough, say ||2;|| > A. Obviously in this 
case, fj5.2j) is nonnegative. 

Now suppose that < A. If x G 5(0, 1) = aJ then h' = and the required result 
holds. Otherwise we write for a; G aJi \ oJ and xq G du (i.e. in our local coordinates 
|x| > 1 and Ixol = 1) : 

\Dxix)fi{b,x,z) - Dxixo)fi{b,xo,z)\ < C\x - y\ 

for a constant C. But the hypothesis (Hs) for fi writes 

Dx{xo)fi{b,xo,z) > C > 0. 
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Then we distinguish two cases {x near 1 or "far" from 1) and get the nonnegativity 
of ()5.2|) in both situations. This can be done by using similar arguments to those 
displayed at the end of the Third Step in the proof of Proposition 14 .4 . in particular, 
a A large enough is needed. 

Third Step : The first two steps give the existence of processes (X', Z') (with 
e a;) solving equation {M + D)r associated to the drift f and terminal value . 
But /' converges to / uniformly on uJ so, passing through the limit as in Subsection 
14.61 we get a pair of processes (Xj, Zt)o<t<r, with X E ZJ and solving the initial 
equation {M -\-D)t- with drift / and terminal value If^. This completes the proof. □ 

We consider again a stopping time r and the corresponding equation (M + D)r] 
now, we only suppose that r is finite a.s. and verifies the exponential integrability 
condition Examples of such stopping times are exit times of uniformly elliptic 

diffusions from bounded domains in Euclidean spaces. 

In this case, we need to add restrictions on the drift /; indeed, the main thrust 
in the proof of uniqueness and existence is the construction of a submartingale on 
the product manifold {St)t = (exp(At)\E'(Xt))j with /i = if / does not depend 
on z. To extend this approach to a random (non necessarily bounded) interval, 
we have to keep the integrability of Sr- An accurate examination of the method 
to obtain the submartingale (in particular inequalities ()3.5|) and ()3.9|) ) shows that 
this integrability holds for "small" drifts; more precisely there is a constant h with 
< h < p such that, under the following condition on the constants in p.4|) and 



the integrability required holds, so {St)o<t<T is a true submartingale. 

Remarks : 1- Such a condition guarantees in particular that we have 



2- This condition is rather natural; actually, it is very similar to conditions yet 
introduced for BSDEs with Lipschitz coefficients and random terminal time : see 
condition (24) and Propositions 3.2 and 3.3 in T], or (2.6) and the condition before 
in [26J. 

3- A priori, the process {Zt)t verifies the integrability condition 



in fact, it results from the existence part that in any case (i.e. / depending or not 
on z), (^t)o<t<r belongs to (So) (see Definition 13.4. 1|) for a small enough, which is 
a stronger property. In particular, we get 



L < h, L2 < h, 



(5.3) 
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this condition is usual for BSDEs with random terminal time (see again Propositions 
3.2 and 3.3 in |7j, or Theorem 2.2 in ^^). 

Once we have constructed the submartingale as on a deterministic interval, 
uniqueness is straightforward. Let us indicate how existence can be deduced. 

We are given a uJ- valued and JF^- measurable variable U"^ . As in the proof of 
Theorem 15.3. !( we consider again local coordinates introduced in Subsection 14.41 ZJi 
such that tJ C tJi C O and a cut-off function with = 1 on aJ and = outside 
ZJi. We put X, z) = (pix)f{b, x, z). 

The first step here is to solve on [0; r] a BSDE whose terminal value is near : 
on [0; r A n], using Theorem 15.3.11 we solve equation (M + U) with drift /i and 
terminal value at time t f\n, E[f/'^|jF„]; let [X"^ , Z")o<t<T-An denote the solution; 
on [r A n; r], we put Zt = and solve 

dX^ = MBlX^,0)dt 

Then it is easily seen, since E[?7'^|jF„] is jF^An-measurable, that {X^,Z^)Q<t<T is a 
solution to BSDE (M + D)r with terminal value C/"^'", where 



JrAn 



' tAti 

The second step is to show that when n tends to infinity, we get the solution of 

BSDE (M + D)^ with terminal value W. 

We have that U^'^ tends to U"^ in indeed. 



E\U^ - f/^'"P = E 1 



- E[f/^|^„,] - / fi{By,x:,o)ds 



rAn 



2' 



and the last expectation tends to zero as n tends to infinity; this is a consequence 
of dominated convergence, using the exponential integrability condition of r. 
Then we apply the results of Subsection 14.11 to the random interval [0;r]; indeed, 
an accurate examination shows that these results rely essentially on the exponential 
integrability condition (according to the uniqueness part) 



E (^e"/o^(ll^?iP+ll^riP)<^^j < C < oo 



for a constant C independent of m, n. At the end, we get the existence of a pair of 
processes {Xt, Zt)Q<t<T solution of BSDE {M + D)j. with drift /i and terminal value 
U'^ . But, since U'^ is tJ-valued and uJ = B{0, 1), for each n the process (X")o<t<rAn 
remains in uJ by Theorem 15.3.11 thus the whole process {Xt)o<t<T remains in uJ and 
this completes the existence part. 

As a consequence, we can state the following result 

Theorem 5.3.2 We consider BSDE {M + D)r with t a stopping time verifying the 
integrability condition ()l.(i|l .- the function x used to define the domain uJ is supposed 
as usual to be strictly convex. Then if f verifies conditions ()1.4j) . ()1.5j) . (H) and 
moreover is "small" (i.e. verifies condition ()5.3|) above), this BSDE has a unique 
solution {X,Z), in the same cases as in Theorem \1.4.1\ 
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5.4 Application to nonlinear elliptic PDEs 

In this paragraph, we make precise the Dirichlet problem that we briefly discussed 
in the introduction. 

Suppose (A^, g) is a Riemannian manifold, and a Brownian motion on (A^, g) 
(started at x at time 0). Alternatively, think of as the diffusion process on Mf^, 
defined by p.2|) : in this case, 

d 

Vz, J = 1, . . . , d, ^ auCTji = g'^ , 
1=1 

the inverse metric tensor, and 

Wz = l,...,d, + J2 g'^Tl, = 0. 

k,l=l 

Let Ml be a compact sub manifold of A^, with boundary dMi and interior Mi. For 
X e Ml, we call ( the first time B^ hits the boundary; we assume that ( verifies an 
integrability condition like (jl.fij] . Given a regular mapping 

: dMi ^ cJ C M, 

we wish to find a mapping : Mi ZJ which solves the following Dirichlet problem 

, .(CM(p{x)-f{x,(f){x),y4>{x)(j{x)) = , X e Ml 
^ ' \ = 0(x) , X G dMi 

where Cm(1> is the tension field of the mapping cj) (see f^, or for a probabilistic point 
of view the introduction of [^). 

We recall from the introduction that, in coordinates (x*) on M and {y") on Mi, 
the equation Cm4' = characterizes harmonic mappings, and is written 

V^, AMi(l)' + g''^T],{<p)D^<pWf,<p' = 0. 

Using the same Wiener process W with which we constructed B^, we can solve 
according to Theorem I5.H.2I the BSDE (M + D)^ with terminal value (f){B^). Let 
{Xf , Zf)o<t<( be the unique solution and put 0(x) := Xq. Then under sufficient 
regularity on 0, it is not difficult to verify that is a solution to the Dirichlet 
problem {D). Note that when / = (i.e. in the martingale case), Kendall (|17j) 
has proved regularity results on using almost only probability theory, so that is 
a strong solution of the equation Cm4' = (i-G- ^ harmonic mapping). 

When f{b,x,z) = f{b,x) and is written as f{b,x) = D2G{b,x) (the differential 
of G with respect to the second variable), the elliptic nonlinear PDE in the Dirichlet 
problem (D) is associated with a variational problem; more precisely, solutions of 
this equation are critical points of the functional 

J^{u) = \f \\giaduib)fdvol{b)+ [ G{b,u{b))dvol{b) 

2 J Ml J Ml 

and the elliptic PDE in equation {D) is the Euler-Lagrange equation associated. 
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5.5 Application to nonlinear parabolic PDEs 

We conclude this part by studying the time-dependent equation associated with the 
stationary equation described in the Dirichlet problem (D) above. More precisely, 
we are interested in the following equation, for mappings m : [0; T] x iV — > aJ C M : 

^ = Cmu- f{x,u,Vua) 
u\t=o = F 

where F is sufficiently regular and has range uj. In local coordinates, this equation 
becomes 

f (t,x) = \Cu{t,x) + \Tjk{u{t,x)){[{V,ua){t,x)Y\[{V,ua){t,x)]^) 

- f{x,u{t,x),{V^ua){t,x)) 

u{0,x) = F{x) 

in the case of C being the Laplace-Beltrami operator on A^. This is equation ()4.6|) . As 
a by-product of Section EJ we have the existence and uniqueness of a regular solution 
to this system of quasilinear parabolic PDEs; it is based on the boundedness of VxU, 
proved in Subsection 14.31 
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